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ABSTRACT 

Our Chandra X-ray Observatory archival study of intracluster entropy in a sample of 222 galaxy clusters 
shows that Ha and radio emission from the brightest cluster galaxy are much more pronounced when the 
cluster's core gas entropy is < 30 keV cm 2 . The prevalence of Ha emission below this threshold indicates 
that it marks a dichotomy between clusters that can harbor multiphase gas and star formation in their cores 
and those that cannot. The fact that strong central radio emission also appears below this boundary suggests 
that AGN feedback turns on when the intracluster medium starts to condense, strengthening the case for AGN 
feedback as the mechanism that limits star formation in the Universe's most luminous galaxies. 
Subject headings: conduction -cooling flows -galaxies: active -galaxies: clusters: general -X-rays: galaxies: 
clusters - galaxies: evolution 



1. INTRODUCTION 

In recent years the "cooling flow problem" has been the 
focus of intense scrutiny as the solutions have broad impact 
on ou r theories of galaxy formation (see lPeterson"&~ Fabian 
120061 for a review). Current models predict that the most 
massive galaxies in the Universe - brightest cluster galax- 
ies (BCGs) - should be bluer and more massive than obser- 
vations find, unless AGN feedback intervenes to stop late- 
time star forma tion dBower et al.l 120061 : ICroton et alJ 120061 
ISaro et alJl2006h . X -ray observations of galaxy clusters have 
given this hypothesis considerable traction. From the prop- 
erties of X- r ay cav ities in the intracluster medium (ICM), 
iBirzan et all (120041) concluded that AGN feedback provides 
the nece ssary energy to retard coo ling in the cores of clus- 
ters (see lMcNamara & Nu lsen 2007 for a review). This result 
suggests that, under the right conditions, AGN are capable of 
quenching star formation by heating the surrounding ICM. 

If AGN feedback is indeed responsible for regulating star 
formation in cluster cores, then the radio and star-forming 
properties of galaxy clusters should be related to the distribu- 
tion of ICM specific entropy 1 . In previous observational w ork 
(see lDonahue et alJl2005Ll2006l ICavagnolo et alj|2008ah . we 
have focused on ICM entropy as a means for understanding 
the cooling and heating processes in clusters because it is a 
more fundame ntal property of the ICM than temperature or 
density alone (IVoit et alJl2002h IVoitll2005l) . ICM temperature 
mainly reflects the depth and shape of the dark matter poten- 
tial well, while entropy depends more directly on the history 
of heating and cooling within the cluster and determines the 
density distribution of gas within that potential. 

We have therefore undertaken a large Chandra archival 
project to study how the entropy st ructure of clusters corre- 
lates with other cluster properties. Cavagnolo et al. ( 2008a) 
presents the radial entropy profiles we have measured for a 
sample of 222 clusters taken from the Chandra Data Archive. 
We have named this project the Archive of Chandra Cluster 
Entropy Profile Tables, or ACCEPT for short. To characterize 
the ICM entropy distributions of the clusters, we fit the equa- 
tion K(r) = Xo+Xioo(r/100 kpc) Q to each entropy profile. In 
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In this paper we quantify entropy in terms of the adiabatic constant K - 
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this equation, K\oq is the normalization of the power-law com- 
ponent at 100 kpc and we refer to Kq as the central entropy. 
Bear in mind, however, that Kq is not necessarily the minimum 
core entropy or the entropy at r = 0, nor is it the gas entropy 
which would be measured immediately around the AGN or 
in a BCG X-ray coronae. Instead, Kq represents the typical 
excess of core entropy above the best fi tting power-law found 
at larger radii. ICavagnolo et all d2008al) shows that Kq is non- 
zero for almost all clusters in our sample. 

In this letter we present the results of exploring the relation- 
ship between the expected by-products of cooling, e.g. Ha 
emission, star formation, and AGN activity, and the Kq val- 
ues of clusters in our survey. To determine the activity level 
of feedback in cluster cores, we selected two readily avail- 
able observables: Ha and radio emission. We have found 
that there is a critical entropy level below which Ha and radio 
emission are often present, while above this threshold these 
emission sources are much fainter and in most cases unde- 
tected. Our results suggest that the formation of thermal in- 
stabilities in the ICM and initiation of processes such as star 
formation and AGN activity are closely connected to core en- 
tropy, and we suspect that the sharp entro py threshold we 
have found arises from thermal conduction (I Voit et al.ll2.Q08l 
in press). 

This letter proceeds in the following manner: In $2] we 
cover the basics of our data analysis. The entropy-Ha rela- 
tionship is discussed in $3] while the entropy-radio relation- 
ship is discussed in S A brief summary is provided in <J3] 
For this letter we have assumed a flat ACDM Universe with 
cosmogony Qm = 0.3, il\ = 0.7, and Hq = 70 km s" 1 Mpc -1 . 
All uncertainties are 90% confidence. 

2. data analysis 

This section briefly describes our data reduction and meth- 
ods for producing e ntropy profiles . Mor e thorough expla- 
nations are given in Donahue et al] (120061) . ICavagnolo et all 
(2008ah. and TCavagnolo et ail (l2008bl) . 



2.1. X-ray 

X-ray data was taken from publicly available observations 
in the Chandra Data Archive. Following standard CIAO re- 
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duction techniques 2 , data was reprocessed using CIAO 3.4.1 
and CALDB 3.4.0, resulting in point source and flare clean 
events files at level-2. Entropy profiles were derived from the 
radial ICM temperature and electron density profiles. 

Radial temperature profiles were created by dividing each 
cluster into concentric annuli with the requirement of at least 
three annuli containing a minimum of 2500 counts each. 
Source spectra were extracted from these annuli, while cor- 
responding background spectra were extracted from blank- 
sky backgrounds tailored to match each observation. Each 
blank-sky background was corrected to account for variation 
of the hard-particle background, while spatial variation of 
the soft-galactic background was accounted for through ad- 
dition of a fixed background component during spectral fit- 
ting. Weighted responses which account for spatial variations 
of the CCD calibration were also created for each observa- 
tion. Spectra were th en fit over the energy range 0.7-7.0 keV 
in XSPEC 1 1.3.2ag ( lArnaudll9 96) using a single-component 
absorbed thermal model. 

Radial electron density profiles were created using surface 
brightness profiles and spectroscopic information. Exposure- 
corrected, background-subtracted, point-source-clean surface 
brightness profiles were extracted from 5" concentric annular 
bins over the energy range 0.7-2.0 keV. In conjunction with 
the spectroscopic normalization and 0.7-2.0 keV count rate, 
surface brightness was co nverted to electron density using the 
deprojection technique of Rri ss et all (119831) . Errors were es- 
timated using 5000 Monte Carlo realizations of the surface 
brightness profile. 

A radial entropy profile for each cluster was then produced 
from the temperature and electron density profiles. The en- 
tropy profiles were fit with a simple model which is a power- 
law at large radii and approaches a constant value, Kq, at small 
radii (see |l]for the equation). We define central entropy as 
Ko from the best-fit model. 

2.2. Ha 

One goal of our project was to determine if ICM en- 
tropy is connected to processes like star formation. Here 
we do not directly measure star formation but instead use 
Ha, which is usuall y a strong indicat or of ongoing star for- 
mation in galaxies (Kennicutil 119831) . It is possible that 
some of the Ha emission from BCGs is not produced by 
star formation (j Begelman & Fabian 1990; Sparks et~aT]l2004t 
iRuszkowski et all 120081: iFerland et alj 120081) . Nevertheless. 
Ha emission unambiguously indicates the presence of ~ 10 4 
K gas in the cluster core and therefore the presence of a mul- 
tiphase intracluster medium that could potentially form stars. 

Our Ha val ues have been gathered from several sources, 
most notably Craw ford et alJ (Il999h . Additional sources 
of data are M. Donahue's observation s taken at Las 
Campa na s and Paloma r (rep ort ed in ICavagnolo et al.l 
2008al). iHeckman et alJ (1 1 9891). iDonahue et al I d!992h . 
Lawrence et al J dl9%h. IValluri & Anupamal (1 19961) 
White et all dl997l) . ICrawford et al.l d2005l) . and lOuillen et alJ 
(2008). We have recalculated the Ha luminosities from these 
sources using our assumed ACDM cosmological model. 
However, the observations were made with a variety of 
apertures and in many cases may not reflect the full Ha 
luminosity of the BCG. The exact levels of Lu a are not 
important for the purposes of this letter and we use the Lhq 
values here as a binary indicator of multiphase gas: either Ha 

2 http://cxc.harvard.edu/ciao/guides/ 



emission and cool gas are present or they are not. 

2.3. Radio 

Another goal of this work was to explore the relation- 
ship between ICM entropy and AGN activity. It has long 
been known that BCGs are more like ly to host radio- 
loud AGN than other cluster galaxies dBurns et al.l 119811 : 
IValentijn & BijleveldlfT983l:lBurnslll990l) . Thus, we chose to 
interpret radio emission from the BCG of each ACCEPT clus- 
ter as a sign of AGN activity. 

To make the radio measurements, we have taken advan- 
tage of the nearly all-sky fl ux-limited coverage of the NRAO 
VLA Sky Survey (NVSS, ICondon et all [1991 and Sydney 
University Molong lo Sky Survey (SUMSS lBock et al.lll999l: 
iMauch et all 120031) . NVSS is a continuum survey at 1.4 
GHz of the entire sky north of S = -40°, while SUMSS is 
a continuum survey at 843 MHz of the entire sky south of 
S = -30°. The completeness limit of NVSS is « 2.5 mJy and 
for SUMSS it is m 10 mJy when 6 > -50° or « 6 mJy when 
5 < -50°. The NVSS positional uncertainty for both right 
ascension and declination is < 1" for sources brighter than 
15 mJ y, and rj 7" at the survey detection limit (Cond on" et alJ 
1998). At z = 0.2, these uncertainties represent distances on 
the sky of ~ 3-20 kpc. For SUMSS, the positional uncer- 
tainty is < 2" for sources brighter than 20 mJy, an d is always 
less than 10" dBock et al.lU999t IMauch et al.ll2003l) . The dis- 
tance at z = 0.2 associated with these uncertainties is ^6-30 
kpc. We calculate the radio power for each radio source us- 
ing the standard relation vh v = AiiD\S u fQ where S v is the 1 .4 
GHz or 843 MHz flux from NVSS or SUMSS, D L is the lu- 
minosity distance, and /o is the central beam frequency of the 
observations. Our calculated radio powers are simply an ap- 
proximation of the bolometric radio luminosity. 

Radio sources were found using two methods. The first 
method was to search for sources within a fixed angular dis- 
tance of 20" around the cluster X-ray peak. The probability of 
randomly finding a radio source within an aperture of 20" is 
exceedingly low (< 0.004 for NVSS). Thus, in 222 total field 
searches, we expect to find no more than one spurious source. 
The second method involved searching for sources within 20 
projected kpc of the cluster X-ray peak. At z ~ 0.051, 1" 
equals 1 kpc, thus for clusters at z > 0.05, the 20 kpc aper- 
ture is smaller than the 20" aperture, and the likelihood of 
finding a spurious source gets smaller. Both methods produce 
nearly identical lists of radio sources with the differences aris- 
ing from the very large, extended lobes of low-redshift radio 
sources such as Hydra A. 

To make a spatial and morphological assessment of the ra- 
dio emission's origins, i.e. determining if the radio emission 
is associated with the BCG, high angular resolution is neces- 
sary. However, NVSS and SUMSS are low-resolution surveys 
with FWHM of « 45". We therefore cannot distinguish be- 
tween ghost cavities/relics, extended lobes, point sources, re- 
accelerated regions, or if the emission is coming from a galaxy 
very near the BCG. We have handled this complication by vi- 
sually inspecting each radio source in relation to the optical 
(using DSS I/II 3 ) and infrared (using 2MASS 4 ) emission of 
the BCG. We have used this method to establish that the radio 
emission is most likely coming from the BCG. When avail- 
able, high resolution data from VLA FIRST 5 was added to the 

3 http://archive.stsci.edu/dss/ 

4 http://www.ipac.caltech.edu/2mass/ 

5 http://sundog.stsci.edu 
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FIG. 1 . — Central entropy vs. Ha luminosity. Orange circles represent Ha 
detections, black circles are non-detection upper limits, and blue squares with 
inset red stars or orange circles are peculiar clusters which do not adhere to 
the observed trend (see text). A215 1 is plotted using the 2cr upper-limit of the 
best-fit Kq and is denoted by a green triangle. The vertical dashed line marks 
Kq = 30 keV cm 2 . Note the presence of a sharp Ha detection dichotomy 
beginning at Kq < 30 keV cm 2 . 



visual inspection. VLA FIRST is a 10,000D° high-resolution 
(5") survey at 20 c m of the North and South Galactic caps 
( Beck eret al.ll 1 9951) . FIRST is also more sensitive than either 
NVSS or SUMSS with a detection threshold of 1 mJy. 

3. Ha EMISSION AND CENTRAL ENTROPY 

Of the 222 clusters in ACCEPT, we located Ha observa- 
tions from the literature for 110 clusters. Of those 110, Ha 
was detected in 46, while the remaining 64 have upper lim- 
its. The mean central entropy for clusters with detections is 
Kq =13.9 ±4.9 keV cm 2 , and for clusters with only upper- 
limits K = 130 ± 55 keV cm 2 . 

In Figure Q] central entropy is plotted versus Ha luminos- 
ity. One can immediately see the dichotomy between clusters 
with and without Ha emission. If a cluster has a central en- 
tropy < 30 keV cm 2 then Ha emission is usually "on", while 
above this threshold the emission is predominantly "off". For 
brevity we refer to this threshold as ^thresh hereafter. The clus- 
ter above ^thresh which has Ha emission (blue square with in- 
set orange circle) is Zwicky 2701 (Kq = 39.7 ±3.9 keV cm 2 ). 
There are also clusters below AWesh without Ha emission 
(blue squares with red stars): A2029, A2107, EXO 0422-086, 
and RBS 533. A2151 also lies below /^thresh an d nas no de- 
tected Ha emission, but the best-fit Ko for A2151 is statis- 
tically consistent with zero and this cluster is plotted using 
the 2a upper-limit of Kq (green triangle in Fig. [TJ. These 
five clusters are clearly exceptions to the much larger trend. 
The mean and dispersion of the redshifts for clusters with and 
without Ha are not significantly different, z = 0.124 ±0.106 
and z = 0.1 32 ±0.084 respectively, and applying aredshift cut 
(i.e. z = 0-0.15 or z = 0.15-0.3) does not change the Kq- 
Ha dichotomy. Most important to note is that changes in the 
Ha luminosities because of aperture effects will move points 
up or down in Figure Q] while mobility along the Kq axis is 
minimal. Qualitatively, the correlation between low central 
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FIG. 2. — BCG radio power vs. Kq for clusters with z < 0.2. Orange sym- 
bols represent radio detections and black symbols are non-detection upper- 
limits. Circles are for NVSS observations and squares are for SUMSS obser- 
vations. The blue squares with inset red stars or orange circles are peculiar 
clusters which do not adhere to the observed trend (see text). Green triangles 
denote clusters plotted using the 2<r upper- limit of the best-fit Kq . The vertical 
dashed line marks Ko = 30 keV cm 2 . The radio sources show the same trend 
as Ha: bright radio emission is preferentially "on" for J$ 30 keV cm 2 . 



entropy and presence of Ha emission is very robust. 

The clusters with Ha detections are typically between 
10-30keV cm 2 , have short central cooling times (< 1 Gyr), 
and under older nomenclature would be classified as "cool- 
ing flow" clusters. It has long been known that star for- 
mation and associated Ha nebulosity appear only i n cluster 
cores with cooling times le s s than a Hubble time dHu et alj 

19891: 



1981 Johnston e eU^ [1987: McNamar a~& O' Connelll 1 1 98S 
Voit & Do nahue 1997; Cardiel et al. 1991). However, our re- 
sults suggest that the central cooling time must be at least a 
factor of 10 smaller than a Hubble time for these manifesta- 
tions of cooling and star formation to appear. It is also very 
interesting that the charac teristic en t ropy t hreshold for strong 
Ha emission is so sharp. IVoit et all d2008l) have recently pro- 
posed electron thermal conduction may be responsible for set- 
ting this threshold. This hypothesi s has rece i ved f urther sup- 
port from the theoretical work of iGuo e t al. (2008) showing 
that thermal conduction can stabilize non-cool core clusters 
against the formation of thermal instabilities, and that AGN 
feedback may be required to limit star formation when con- 
duction is insufficient. 

4. RADIO SOURCES AND CENTRAL ENTROPY 

Of the 222 clusters in ACCEPT, 100 have radio-source de- 
tections with a mean Kq of 23.3 ± 9.4 keV cm 2 , while the 
other 122 clusters with only upper limits have a mean Kq of 
134±52keV cm 2 . NVSS and SUMSS are low resolution sur- 
veys with FWHM at « 45" which at z = 0.2 is w 150 kpc. 
This scale is larger than the size of a typical cluster cooling 
region and makes it difficult to determine absolutely that the 
radio emission is associated with the BCG. We therefore fo- 
cus only on clusters at z < 0.2. After the redshift cut, 135 
clusters remain - 64 with radio detections (mean Kq= 18.3 ± 
7.7 keV cm 2 ) and 71 without (mean K = 1 12 ±45 keV cm 2 ). 
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In Figure[2]we have plotted radio power versus Kq. The ob- 
vious dichotomy seen in the Ha measures and characterized 
by Xthresh> is a ls° present in the radio. Clusters with vL v > 
10 40 erg s " 1 generally have Kg < .Kthresh- This trend was first 
evident in iDonahue et al.l (l2005h and suggests that AGN ac- 
tivity in BCGs, while not exclusively limited to clusters with 
low core entropy, is much more likely to be found in clusters 
which have a core entropy less than K± KS h. That star forma- 
tion and AGN activity are subject to the same entropy thresh- 
old suggests the mechanism which promotes or initiates one 
is also involved in the activation of the other. If the entropy 
of the hot gas in the vicinity of the AGN is correlated with 
Ko, then the lack of correlation between radio power and Ko 
below t he 30 keV cm 2 threshold suggests that cold-mo de ac- 
cretion (Pizzolato & Soker 2005; Hardcastle et al. 2007) may 
be the dominant method of fueling AGN in BCGs. 

We have again highlighted exceptions to the general trend 
seen in Figure [2] clusters below /T t hi e sh without a radio 
source (blue squares with inset red stars) and clusters above 
^thresh with a radio source (blue squares with inset orange cir- 
cles). The peculiar clusters below Xthresh are A133, A539, 
A 1204, A2107, A2556, AWM7, ESO 5520200, MKW4, MS 
J0440.5+0204, and MS Jl 157.3+5531. The peculiar clusters 
above are 2PIGG J001 1.5-2850, A193, A586, A2063, 
A2147, A2244, A3558, A4038, and RBS 461. In addition, 
there are three clusters, A2151, AS405, MS 0116.3-0115, 
which have best-fit Ko statistically consistent with zero and 
are plotted in Fig. fusing the 2a upper-limit of Ko (green 
triangles). All three of these clusters have detected radio 
sources. 

Finding a few clusters in our sample without radio sources 
where we expect to find them is not surprising given that AGN 
feedback could be episodic. However, the clusters above 
^thresh with a central radio source are interesting, and may 
be special cases of BCGs with embedded coronae. ISun et all 



d2007l) extensively studied coronae and found they are like 
"mini-cooling cores" with low temperatures and high densi- 
ties. Coronae are a low-entropy environment isolated from 
the high-entropy ICM and may provide the conditions neces- 
sary for gas cooling to proceed. And indeed, 2PIGG 0011, 
A193, A2151, A2244, A3558, A4038, and RBS 461 show 
indications that a very c ompact (r < 5 kpc) X-ra y source is 
associated with the BCG dCavagnolo et al.ll2008al) . 

5. SUMMARY 

We have presented a comparison of ICM central entropy 
values and measures of BCG Ha and radio emission for 
a Chandra archival sample of galaxy clusters. We find 
that below a characteristic central entropy threshold of 
Ko ~ 30 keV cm 2 , Ha and bright radio emission are more 
likely to be detected, while above this threshold Ha is not 
detected and radio emission, if detected at all, is significantly 
fainter. The mean Kq for clusters with and without Ha detec- 
tions are K = 13.9 ±4.9 keV cm 2 and#o= 130±55 keV cm 2 , 
respectively. For clusters at z < 0.2 with BCG radio emission 
the mean K = 18.3 ± 7.7 keV cm 2 , while for BCGs with 
only upper limits, the mean Ko = 112 ±45 keV cm 2 . While 
other mechanisms can produce Ha or radio emission besides 
star formation and AGN, if one assumes the Ha and radio 
emission are coming from these two feedback sources, then 
our results suggest the development of multiphase gas in 
cluster cores (which can fuel both star formation and AGN) 
is strongly coupled to ICM entropy. 
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